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This paper discusses the cumulative probability distributions of field strength for 
four 200 km VHF paths in Illinois in terms of a classification of refractive index profiles. 
It is shown that extended elevated layers produce signal enhancements of 10 to 25 db above 
the level observed in unstratified conditions. Assuming the layer characteristics given 
by radiosonde data, the best agreement between calculated and measured values of field 
strength is obtained using a layer model with a linear n-profile. 

The possible influence of smaller layers is also discussed in relation to the observed 
results for conditions judged to be unstratified or well mixed on the basis of sonde data. 



1. Introduction 

It is evident from the discussion in part I that 
the further development of radiometeorological 
parameters would be assisted by a better under- 
standing of the propagation mechanism on non- 
optical paths. In particular, the influence of 
thermal stability on signal level, fading rale and 
wavelength dependence is an important topic re- 
quiring further study. 

The effect of varying meteorological conditions 
on signal characteristics is especially marked in the 
case of paths of "intermediate" length. On such a 
path the radio field strength in the absence of stable 
layers or surface ducts will often be comparable 
with the predicted diffracted field. Furthermore, at 
frequencies up to say 300 Mc/s (X>1 m) relatively 
strong fields will frequently be observed in con- 
ditions favorable to the production of temperature 
inversions in the first 2 km or so above the earth's 
surface. It is the purpose of this paper to discuss 
some aspects of radiometeorology relevant to this 
situation, especially the field strength distribution 
observed on a 200 km path at frequencies between 



72 and 180 Mc/s. The signal characteristics are 
analyzed in terms of a classification of refractive 
index profiles, with the objective of clarifying the 
relative importance of different propagation mech- 
anisms and their influence on the measured field 
st rength dist ribution. '• ,1 

'Fable 1 lists characteristic profile types, the 
assumed mechanism associated with each type, and 
typical meteorological conditions. Selected refer- 
ences are given for each category, and special 
mention may be made here of recent work [du Castel, 
Misme, and Voge, 1960] in France which, to a large 
extent, unifies and extends earlier analyses based 
on the separate concepts of "reflection" and 
"scattering." 

2. Radio and Meteorological Data 

The analysis to follow is limited to paths between 
Chicago and Urbana, 111. (fig. 1), since several 
years of radio data are available for four separate 
wavelengths between 1.67 m and 4. IS m. More- 
over, two radiosonde stations are located on or near 
the path, a unique situation in radiometeorological 



Table 1. — Refractive index profile classification, probable propagation mechanisms and meteorological conditions 



Profile 



Assumed propagation 
mechanism 



Reference 



Meteorological conditions 



Unstratified, (U) 
Monotonic decrease with height, 
gradient nowhere exceeds twice nor- 
mal for that height. 

Elevated Layer, (EL) 
Monotonic decrease with height with 
one or more distinct layers with gradi- 
ents at least twice normal for that 
height. 

Super-Refractive, (SR) 
Same as EL but the layer is ground- 
based. 

Ducting, (D) 
Same as SR but the gradient exceeds 
the earth's curvature, 1/a. 



Scattering plus diffraction. 



Scattering plus 
Diffraction plus reflection. 



Extended radio horizon pro- 
ducing enhanced diffracted 
and scattered components. 

Extension of radio horizon to 
include the receiver. 



[Booker and Gordon, 1951] 

[Villars and Weiskopf, 1954, 

1955]. 
[Norton et al, 1955]. 

[Smyth and Trolese, 1948] 

[Saxton, 1951], 

[Misme, duCastel, Voge, 19G0] 



[Norton et al, 1955].. 



[Booker- Walkinshaw, 1945].... 



Well-mixed atmosphere due 
to thermal convection and, 
or wind shear. 

Layer formed by subsidence 
inversion or lifting of radia- 
tion inversion. 



Radiation inversion formed 
during the night or rapid 
evaporation from soil after 
rain. 

Same as SR. 



597 




MISSOURI *><""*"- 



Figure 1. Location of radio path and radiosonde stations 
used in this study. 

investigations. Details of the radio paths are 
given in table 2 (in which 6 is the total angle be- 
tween the horizon rays from transmitter and receiver 
on a 4/3-earth profile). 

The meteorological data were obtained from the 
simultaneous radiosonde observations made by the 
Weather Bureau at Joliet and the Rantoul Air 
Force Base. The results used were those from the 
significant levels reported whenever the temperature 
or humidity departed by 1°C or ±10 percent from 
predetermined values. 



Table 2- 


—Chicago-to-Urbana radio path characte? istics 


Station 


Distance 





X 


/ 


Period of 
record 


WBKB-TV 


km 

203.1 

202.9 

202.7 

202.9 


mr 
16.3 
16.7 
16.1 

16.4 


m 

4.18 
3.67 
3.15 

1.67 


Mc/s 
71.75 
81.75 
95.50 

179. 75 


5/51- 5/53 


WNBQ-TV 


10/50-10/51 


WMBI-FM 


7/50- 6/52 


WENR-TV 


7/51- 6/53 







2.1. Classification of Radio Field Strengths by Profile 
Types 

The RAOB significant level data were converted 
to refractive index by use of the Smith-Weintraub 
[1953] relationship 



JV= (72,-1) 10 



T V 



P- 



4810A 
T ) 



(1) 



where the pressure, P, temperature, T, and vapor 
pressure, e, are in the usual units of mb, °K, and mb 
respectively. The gradient of N was then de- 
termined between the reported significant levels of 
each profile and examined as to whether the gra- 
dients fell into the category of linear, subrefractive 
or superrefractive depending upon the criteria set 
down in table 3 wherein superrefractive is approxi- 
mately twice normal and subrefractive has a positive 
gradient. Simultaneous observations of similar pro- 
file types at Rantoul and Joliet were necessary for 
entry as a distinct profile occurrence. If a super- 
refractive layer occurred above the crossover heights 
of the radio horizon tangent rays from both trans- 
mitter and receiver, then it was classified as an ele- 
vated layer provided the reported layer heights 
were within 1 km of one another at both radiosonde 
stations. Elevated layers below the crossover height 
were classified as ground-based superrefractive layers. 
Elevated layers below the crossover height at one 
weather station and above that height at the other 
were classified as tilted elevated layers. 

After these characteristic profiles were isolated, 
the median field strengths for the 3 hr period centered 
upon the radiosonde observation time were arranged 
into cumulative probability distributions for each 
profile type. The results are shown in figure 2. 
(There were relatively few examples of subrefractive 
profiles and no distributions for this category are 
given.) The 3-hr time interval about the 10 am 
and 10 pm CST radiosonde observation times was 
arbitrarily chosen to smooth the sometimes abrupt 
changes in signal characteristics observed during 
these transition periods of the diurnal cycle. 



Table 3. N gradient classification of profile types in N- 
units/km 



p 


ft 


Sub- 
refractive 


Unstratifled 


Super- 
refractive 


mb 
1000-850 

850-700 

700-600 

600-500 

500-400 


km 
-1.46 

1.46-3.01 

3. 01-4. 20 

4. 20-5. 57 

5. 57-7. 18 




20<-^<60 

an 

20<-^<50 
an 

15<-^<40 

an 

10 <_^<30 

dh 

10<-^< 25 
an 


»<-S 

"Hi 
•<-S 



Generally, the unstratifled samples have the lowest 
overall field strengths throughout the entire distri- 
bution range. The presence of any layer (elevated 
or ground-based) tends to increase the field strength 
by 10 to 25 db at any percentage level of the distri- 
bution. (The exception to this observation, WNBQ- 
TV, is probably explained by the fact that the obser- 
vations were limited to six winter months, rather 
than the 2-yr period of the other stations.) Tilted 
elevated layers appear to produce the greatest 
enhancement of signal strength, probably as a 
result of focusing effects due to the layer tilt. 
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Figure 2. 



Distribution of hour!// median field strengths with different radiosonde profile conditions. 



2.2. Prediction of Field Strength for Unstratified 
Conditions 

The field strengths recorded during the times 
when the radiosonde ascents at both Joliet and 
Rantoul indicated nonstratified conditions were 
compared with the values predicted by Norton, 
Rice, and Vogler [1955] for the case of diffracted 
plus scattered radio waves. This particular predic- 
tion process is adjustable for the average refractive 
conditions over the path in that it adjusts the 
effective earth's radius factor to the initial gradient 
of N for the calculation of diffracted field strengths. 
One also needs the angular separation of the radio 



horizon rays at their intersection near midpath. 
The average initial gradient of N was obtained for 
each instance of unstratified profile by simply aver- 
aging the initial gradients from Rantoul and Joliet, 
while the angular separation was obtained by deter- 
mining the amount of radio ray refraction expected 
over each particular path in atmospheres of expo- 
nential decrease with height that closely match the 
observed N conditions. 

Figure 3 illustrates a comparison of the predicted 
and observed field strengths. For WNBQ and 
WENR there is approximate agreement between 
the two sets of data. However, the predicted values 
for WBKB and WMBI are approximately 10 to 12 
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Figure 3. Comparison of observed radio field strengths and values predicted for wintertime afternoon hours. 



db higher than the observed values. This tendency 
to predict fields in excess of the measured values 
suggests that the empirical data, on which the 
predictions are based, may include meteorological 
conditions with some degree of stratification, even 
though most of the empirical data refer to afternoon 
hours in winter. It will be shown in section 3.2 
that elevated layers of moderate size (say a few 
kilometers in horizontal dimensions) may exist 
undetected by the radiosonde and could produce field 
strengths on the Illinois paths comparable with the 
median values for "unstratified" conditions shown 
in figure 2. Furthermore, the limits placed on the 
profile gradients specifying unstratified conditions in 
table 3 are such that some layer type profiles may 
be included in the unstratified category. Conse- 
quently, it is important to study in more detail the 
properties of the elevated layer, not only as a feature 
occasionally producing high field strengths, but also 
as a mechanism which, in less intense form, partly 
determines the strength of the weaker fields observed 
for large percentages of the time. 



3. Effect of Elevated Layers on the Illinois 
Paths 

The influence of elevated layers on VHF trans- 
mission beyond the horizon has been studied by 
several workers [Saxton, 1951; Gossard and Ander- 
son, 1956; du Castel, Misme, and Voge, 1960; Friis, 
Criawford, and Hogg, 1957; Starkey et al., 1958 
Abld et al., 1952]. However, few investigations; 
have contained any detailed comparisons of theory 
and experiment results. The following analysis pre- 
sents such a comparison, using simple models of the 
elevated layer, for the four Illinois paths. 

3.1. Elevated Layers at Temperature Inversions 

Recent radar and refractometer investigations of 
tropospheric structure have shown that elevated 
layers in the refractive index distribution are fre- 
quently observed in the stable air of temperature 
inversions [Lane and Meadows, 1963]. A typical 
value of layer thickness is 100 m, with horizontal 
dimensions of tens of kilometers. On occasions, 
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extended layers no more than 10 m in thickness have 
been detected by refractometer soundings. In the 
present discussion we attempt to evaluate the re- 
flection coefficient of these elevated layers. We 
may express the modulus of the reflection coefficient 
|p|, for a wave incident at a glancing angle a on a 
layer of thickness h, in the form: 



=£f(a,h,z). 



(1) 



f(a, h, x) is the ratio of the reflection coefficient of 
the model to that of the infinitely sharp case (i.e., 
the Fresnel discontinuity value, An/aa 2 ). This 
function has been evaluated for several layer profiles 
[du Castel, Misme, and Voge, 1960], and preliminary 
calculations based on this work were made to de- 
termine the most suitable model in the present 
application. It was evident from these calculations 
that a simple linear profile would yield the best 
agreement with the measured data, and this model 
was therefore adopted in the subsequent analysis. 
Consider the layer profile shown in figure 4, i.e., 
a linear decrease of n over a height interval h, with 
transition regions of height d. This model and others 
have been discussed by several authors, but the 
most detailed treatment is that of Brekhovskikh 
[I960.] His analysis shows that for this linear 
model: 

\p\ = An'\I^Trh sin 3 a 



c^An-X/Sirha 2 . 



(2) 



X 




FlGUKE 4. 



REFRACTIVE INDEX 

Linear profile model of the n decrease across an 
elevated layer. 



This equation is valid if: 

(a) A?i'\<<7rAa 3 
(b) 4o*Z«X. 



and 



In the present problem, with values of A of 1.7 to 
4.2 in, a ~ 10~ 5 , condition (a) is satisfied for layer 
thicknesses greater than about 20 m. In addition, 
condition (b) is fulfilled for the stated conditions 
if the thickness of the transition region is less than 
a few meters. These conditions do not seem incon- 
sistent with available refractometer data on elevated 
layers, but a rigorous justification of the model is 
impossible at the present time. In any case, there 
is almost certainly no unique profile representative 
of all elevated layers. We assume here, therefore, 
the linear profile of figure 4 merely as a simple 
analytical model. It may be noted here that the value 
of \p\ given by (2) agrees with that quoted by 
du Castel [1961], but is half the value obtained in 
an earlier analysis Idu Castel, Misme, and Voge, 
I960]. 

Equation (2) was used to calculate the reflection 
coefficient of the layers on each occasion on which 
these were observed in the sonde ascents. The 
results, expressed in terms of a reflection loss, are 
compared with the measured values of field strength 
in figure 5. The general agreement is satisfactory 
for the assumed model. As might be expected, there 
is a considerable scatter in the data, and two con- 
siderations are important in assessing the signifi- 
cance of these results. These concern sonde 
response and layer structure. The work of Wagner 
[1960] on the response of radiosondes show that, 
for an elevated inversion layer with Aft=3X 10~ 5 , and 
^=100 m, a sonde with a 10 sec time constant in 
the sensing elements, rising at 5 m/sec, will give an 
indicated value of An of approximately half the true 
value. The above procedure, using sonde data, 
therefore underestimates the value of \p\ for an 
idealized infinite layer. On the other hand the 
analysis assumes a smooth layer extending hori- 
zontally at least over a distance x equal to the first 
Fresnel zone. We have: 



:V2aX/2 



(3) 



where 2a = path length = 2-3 km for the Illinois 
paths. Hence x is of the order of a few tens of 
kilometers. In addition, we have assumed that the 
layer is horizontal and smooth over a distance x, 
thus neglecting convergence. If we adopt the 
Rayleigh criterion, the height of the surface irreg- 
ularities on the layer, Ah say, must not exceed ±X/ 
8a for the layer to be considered smooth; i.e., 

Ah<±7 m (X=1.67 m;/=179.75 Mc/s) 

Ah<±17 m (X=4.18 m;/=71.75 Mc/s). 

These values apply for a =0.03 radian, correspond- 
ing to a layer height of 2.5 km; for lower layers 
Ah will be greater due to the decrease in a. 
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Figure 5. Observed hourly median field strength versus reflection loss assuming the linear model of figure 4- 



These conditions are not likely to be satisfied 
in all the examples studied and the analysis therefore 
overestimates the value of \p\ in this respect. (Some 
discussion of this point has been given by Bauer 
and Meyer [1958].) The limitations of sonde sound- 
ings, and the effects of layer tilt and surface ir- 
regularities therefore provide a partial explanation 
of the scatter of points in figure 5. Further detailed 
measurements of layer structure are obviously 
desirable. 

3.2. Influence of Small Layers 

The above discussion has dealt with the particular 
case of extended elevated layers such as are often 
associated with temperature inversions. However, 
it seems quite possible that these layers are merely 
the more extreme examples of anisotropic irregu- 



larities which are thought to be prevalent in the 
troposphere. There is already some preliminary 
evidence supporting this concept in the results of 
refractometer and radar soundings [Lane and 
Meadows, 1963; Saxton, 1960], and recent theoreti- 
cal work [du Castel, Misme, and Voge, 1960; Friis, 
Crawford, and Hogg, 1957; du Castel, 1961] lias 
developed this approach in detail. The relationship 
of this work to earlier analyses in terms of a "scat- 
tering" model is discussed in the references quoted 
and need not concern us here. For our purpose it is 
sufficient to utilize the essential features of the 
argument as the basis for a simple calculation. 

It seems reasonable to assume that even in an 
atmosphere which sonde data would lead us to 
classify as "well mixed" there are often layer-type 
irregularities. Detailed evidence on the spatial 
form and stability of this type of layer or "feuillet" 
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is so far lacking, but an inspection of some refrac- 
tometer results suggests that horizontal dimensions 
of a few kilometers represent a realistic assumption. 
Such a layer might exist as a separate entity for say 
several minutes (as compared with a period of sev- 
eral hours for the extended layer in a stable 
inversion). 

For the following analysis, let us consider two 
layers of horizontal dimensions, x, of 2 and 10 km 
respectively with the following characteristics: 

A7i=10~ 5 
A=100 m 



= 0.01-0.03 radian (i.e., layer 
km on the Illinois paths). 



height of 0.4-2.5 



For these conditions, the layers correspond to those 
of "intermediate" size in the analysis of Friis, 
Crawford, and Hogg. They are defined by the 
equation: 

V2oX/a>z>V2o\ (4) 

where 2a is the path length. In this case, the power 
received, P R , from an antenna of effective aperture, 
A R , with a transmitter radiating a power P T from an 
antenna of effective aperture A T is given by: 



Pr/Pt= (A T A R \ 2 a 2 P 2 )K2\W). 



(5) 



We can use this equation to calculate the correspond- 
ing field strength, for the Illinois paths, in terms of 
/zv/m for 1 kw radiated from a half-wave dipolc. 
We have the following relations: 



A(\/2 dipole)-0.127A 2 
P«(X/2 dipole)=£ T2 X 2 /3007r 2 



(6) 
(7) 



where E is the field strength in volts/meter if P R is 
in watts. From (5), (6), and (7) we can calculate 
E for the two layers specified r^ove, and the results 
obtained are shown in figure 6 for various layer 
heights and the following models of reflection 
coefficient: 

(a) |p|=A7&.\/87ra 3 & 



(b) 



-An/ 2 c 



Model (b) is the Fresnel discontinuity equation 
which gives the limiting value of \p\ toward which 
all models tend as the layer thickness decreases. 
The curves in figure 6 show that the calculated field 
strength depends considerably on the assumed n- 
profile. If \p\ = An-\/STa 3 h, values of field strength 
comparable with the long-term median value may 
be produced by layers of about 10 km in lateral di- 
mensions in the height range 0.5-1 km. If |p| = 
An/2a 2 } similar field strength may be produced by 
layers in this height range if the lateral dimension is 
of the order of 2 km. The effect of the layer de- 
creases with increasing height, but even with layer 
heights of 3 km, the field strength is still 1 pv/m or 
greater at both wavelengths for a 10 km layer with 
\p\ = Anj2a 2 . However, it should be pointed out 
that the assumed value of An=10~ 5 is probably 



somewhat large for layers as high as 3 km. The 
results also show that model (b) (i.e., p=An/2a 2 ) 
gives field strength values which are higher at \= 
1.67 m (f= 179.75 Mc/s) than at X-4.18 m (/= 
71.75 Mc/s). 

The distribution of wavelength dependence ; is 
further illustrated in figure 7. The hourly median 
values of field strength (or transmission loss, L,) at 
times of the sonde ascents were used to derive the 
wavelength dependence as a function of refractive 
index profile characteristics. Here, transmission loss 
is defined by pa/py, where: 

2ta=Available power at the terminals of a loss-free 

antenna. 
py= Power radiated from the transmitting antenna. 
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Figure^ 6. Field 'strength produced by layers of horizontal 
dimensions 2 km and 10 km on Chicago- Urbana paths for 
two model profiles. 
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The wavelength dependence varies considerably, 
even in any one profile category; this result further 
emphasizes the complexity of the propagation mech- 
anism on a VHF path of intermediate length (~200 
km). A full discussion of this problem is outside 
the scope of this paper, but it seems probable that 
all of the following features are important in deter- 
mining the range of wavelength dependence: 

(a) Occasional deep surface or elevated ducts, 
giving enhancements of the 180 Mc/s signal, 

(b) the different profile characteristics of elevated 
layers, as illustrated by the results of figure 6; 

(c) the path geometry, 

(d) elevated layers or large anisotropic eddies 
which are not detected by the radio sondes. 

It seems likely that the weakest fields measured on 
the Illinois paths are produced by a diffracted com- 
ponent. The wavelength dependence in such a sit- 
ation has been discussed by Schelleng, Burrows, and 
Ferrell [1933], who drew attention to enhancements 
of the higher frequency field strength, on a two- 
channel VHF path, resulting from the path geometry 
and consequent antiphase interference between ray 
paths at the lower frequency. Omitting the extreme 
categories of (a) weak diffracted fields, and (b) fields 
influenced by deep ducts, the range of wavelength 
dependence indicates the possible effect of a distri- 
bution of layer sizes with varying surface roughness 
and 7i-profile. Even in an apparently "well-mixed" 
atmosphere there is rarely a unique propagation 
mechanism on this 200 km path in the frequency 
range 70 to 180 Mc/s, a conclusion consistent with 
the results of figure 3. 

4. Conclusions 

Any departure of refractive index structure from 
a smooth monotonic decrease with height produces 
an increase in field strength on a 200-km path in the 
frequency band 70 to 180 Mc/s. In the particular 
case studied, elevated tilted layers result in signal 
enhancements of 10 to 25 db, over the values for 
unstratified conditions, at all percentage levels. 
(The importance of the tilted layer is possibly a 
consequence of the asymmetry of the path, the 
transmitting antenna being 200 m above ground and 
the receiving antenna 30 m.) 

The predicted field strengths, for conditions classi- 
fied as unstratified in terms of sonde data, are in 
approximate agreement with observed results, al- 
though the scatter of the points (plus the tendency 
to predict values in excess of the measured ones) 
point to the influence of anisotropic layers or eddies 
of varying size and degree of stability. This inter- 
pretation is consistent with numerical calculations 
based on the properties of "intermediate" size layers, 



suggested in the analysis of Friis, Crawford, and 
Hogg [1957]. 

Calculations of the field strength produced by 
extended stable layers, using sonde data and a model 
profile with a linear lapse of n with height, are in 
reasonable agreement with the experimental results. 
However, there is probably no unique profile charac- 
teristic of elevated layers. 



The present work has grown out of a study initiated 
some ten years ago by J. W. Herbstreit. His early 
guidance is greatly appreciated. 
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